Measurement of the charged pion electromagnetic form factor 
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Separated longitudinal and transverse structure functions 
for the reaction 1 H(e,e / 7r + )n were measured in the momen- 
tum transfer region Q 2 = 0.6 - 1.6 (GeV/c) 2 at a value of 
the invariant mass W = 1.95 GeV. New values for the pion 
charge form factor were extracted from the longitudinal cross 
section by using a recently developed Regge model. The re- 
sults indicate that the pion form factor in this region is larger 
than previously assumed and is consistent with a monopole 
parameterization fitted to very low Q 2 elastic data. 
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The pion occupies an important place in the study of 
the quark-gluon structure of hadrons. This is exemplified 
by the many calculations that treat the pion as one of 
their prime examples — Q . One of the reasons is that 
the valence structure of the pion, being (qq), is relatively 
simple. Hence it is expected that the value of the four- 
momentum transfer squared Q 2 , down to which a pQCD 
approach to the pion structure can be applied, is lower 
than e.g. for the nucleoli. Furthermore, the asymptotic 
normalization of the pion wave function, in contrast to 



that of the nucleon, is known from the pion decay. 

The charge form factor of the pion, F 7T (Q 2 ), is an essen- 
tial element of the structure of the pion. Rs behaviour at 
very low values of Q 2 , which is determined by the charge 
radius of the pion, has been determined up to Q 2 =0.28 
(GeV/c) 2 from scattering high-energy pions from atomic 
electrons Q . For the determination of the pion form fac- 
tor at higher values of Q 2 one has to use high-energy 
electroproduction of pions on a nucleon, i.e., employ the 
1 H(e,e , 7r + )n reaction. For selected kinematical condi- 
tions this process can be described as quasi-elastic scat- 
tering of the electron from a virtual pion in the proton. 
In the t-pole approximation the longitudinal cross section 
<tl is proportional to the square of the pion form factor. 
In this way the pion form factor has been studied for Q 2 
values from 0.4 to 9.8 (GeV/c) 2 at CEA/Corncll go| and 
for Q 2 = 0.7 (GeV/c) 2 at DESY Q. In the DESY ex- 
periment a longitudinal/transverse (L/T) separation was 
performed by taking data at two values of the electron 
energy. In the experiments done at CEA/Cornell this 
was done in a few cases only, and even then the resulting 
uncertainties in ox were so large that the L/T separated 
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data were not used. Instead for the actual determination 
of the pion form factor ctl was calculated by subtract- 
ing from the measured (differential) cross section a ctt 
that was assumed to be proportional to the total virtual 
photon cross section. No uncertainty in or was included 
in this subtraction. This means that existing values of 
Fit above Q 2 = 0.7 (GeV/c) 2 are not based on L/T sep- 
arated cross sections. This, together with the already 
relatively large statistical (and systematic) uncertainties 
of those data, precludes a meaningful comparison with 
theoretical calculations in that region. 

Because of the excellent properties of the electron 
beam and experimental setup at CEBAF it is now possi- 
ble to determine L/T separated cross sections with high 
accuracy and thus to study the pion form factor in the 
regime of Q 2 = 0.5 - 3.0 (GeV/c) 2 . Using the High Mo- 
mentum Spectrometer and the Short Orbit Spectrometer 
of Hall C and electron energies between 2.4 and 4.0 GeV, 
data for the reaction 1 H(e, e'Tr + )n were taken for cen- 
tral values of Q 2 of 0.6, 0.75, 1.0 and 1.6 (GeV/c) 2 , at a 
central value of the invariant mass W of 1.95 GeV. 

The cross section for pion electroproduction can be 
written as 
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where Ty is the virtual photon flux factor, <f) is the az- 
imuthal angle of the outgoing pion with respect to the 
electron scattering plane and t is the Mandelstam vari- 
able t = (pjr — q) 2 . The two-fold differential cross section 
can be written as 
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where e is the virtual-photon polarization parameter. 
The cross sections ax = ^j? depend on W, Q 2 and t. 
The longitudinal cross section ox is dominated by the t- 
pole term, which contains F v . The 4> acceptance of the 
experiment allowed the interference terms cut and <ttt 
to be determined. Since data were taken at two energies 
at every Q 2 , ol could be separated from or by means of 
a Rosenbluth separation. 

The analysis of the experimental data included the fol- 
lowing p^ |. Electron identification in the Short Orbit 
Spectrometer was done by using the combination of lead 
glass calorimeter and gas Cerenkov containing Freon-12 
at atmospheric pressure. Pion identification in the High 
Momentum Spectrometer was largely done using time 
of flight between two scintillating hodoscope arrays. A 
small contamination by real electron-proton coincidences 
at the highest Q 2 setting was removed by a single beam- 
burst cut on e — 7T+ coincidence time. Then Q 2 , W, 
t, and the mass of the undetected neutron were recon- 



structed. Cuts on the latter excluded additional pion pro- 
duction. Backgrounds from the aluminum target window 
and random coincidences were subtracted. Yields were 
determined after correcting for tracking efficiency, pion 
absorption, local target-density reduction due to beam 
heating, and dead times. Cross sections were obtained 
from the yields using a detailed Monte Carlo (MC) simu- 
lation of the experiment, which included the magnets, 
apertures, detector geometries, realistic wire chamber 
resolutions, multiple scattering in all materials, recon- 
struction matrix elements, pion decay, muon tracking, 
and internal and external radiative processes. 

Calibrations with the overdetermined 1 H(e, e'p) reac- 
tion were critical in several applications. The beam mo- 
mentum and the spectrometer central momenta were de- 
termined absolutely to 0.1%, while the incident beam 
angle and spectrometer central angles were absolutely 
determined to better than 1 mrad. The spectrometer 
acceptances were checked by comparison of data to MC 
simulations. Finally, the overall absolute cross section 
normalization was checked. The calculated yields for e+p 
elastics agreed to better than 2% with predictions based 
on a parameterization of the world data [ fl3| . 

In the pion production reaction the experimental ac- 
ceptances in W, Q 2 and t were correlated. In order to 
minimize errors resulting from averaging the measured 
yields when calculating cross sections at average values 
of W, Q 2 and t, a phenomenological cross section model 
Jl2[ was used in the simulation program. In this cross 
section model the terms representing the ax of Eqn. (|^) 
were optimized in an iterative fitting procedure to glob- 
ally follow the t— and Q 2 -dependence of the data. The 
dependence of the cross section on W was assumed to 
follow the phase space factor (W 2 — M 2 )~ 2 . 

The experimental cross sections can then be calculated 
from the measured and simulated yields via the relation 
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This was done for five bins in t at the four Q 2 -values. 
Here, (Y) indicates that the yields were averaged over the 
W and Q 2 acceptance, W and Q 2 being the acceptance 
weighted average values for that t-bin. Even while the 
average values of W and Q 2 differed slightly at high and 
low e, the use of Eq. (||) with a MC cross section that 
globally reproduces the data allows one to take a common 
average (W, Q 2 ) value. 

A representative example of the cross section as func- 
tion of <j) is given in Figure |l|. The dependence on <j) was 
used to determine the interference terms olt and cttt 
after which the combination cr uns = <7t + was ob- 
tained at both the high and low electron energy in each 
t bin for each Q 2 point. The statistical uncertainty in 
these cross sections ranges from 2 to 5%. Furthermore, 
there is a total systematic uncertainty of about 3%, the 
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FIG. 1. 



most important contributions being: simulation of the 
detection volume (2%), dependence of the extracted cross 
sections on the MC cross section model (typically less 
than 2%), target density reduction (1%), pion absorption 
(1%), pion decay (1%), and the simulation of radiative 
processes (1%) [Q. Since the same acceptances in W 
and Q 2 and the same average values W and Q 2 were 
used at both energies, ox, and <7t could be extracted via 
a Rosenbluth separation. 

These cross sections arc displayed in Figure The 
error bars represent the combined statistical and system- 
atic uncertainties. Since the uncertainties that are un- 
correlated in the measurements at high and low electron 
energies are enlarged by the factor l/(Ae) in the Rosen- 
bluth separation, where Ae is the difference (typically 
0.3) in the photon polarization between the two mea- 
surements, the total error bars on ctl are typically about 
10%. 

The experimental data were compared to the results 
of a Regge model by Vanderhaeghen, Guidal and Laget 
(VGL) 14 1 . In this model the pion electroproduction 
process is described as the exchange of Regge trajecto- 
ries for 7T and p like particles. The only free parameters 
are the pion form factor and the 71737 transition form fac- 
tor. The model globally agrees with existing pion photo- 
and electroproduction data at values of W above 2 GeV. 
The VGL model is compared to the data in Figure ||. The 
value of Fir was adjusted at every Q 2 to reproduce the <tl 
data at the lowest value of t. The transverse cross section 
ctt is underestimated, which can possibly be attributed 
to resonance contributions at W = 1.95 GeV that are 
not included in the Regge model. Varying the irp'f tran- 
sition form factor within reasonable bounds changes ctt 
by up to 30%, but has a negligible influence on <7l, which 
is completely determined by the 7r trajectory. This t- 
pole dominance was checked by studying the reactions 
2 H(e, e'ir + )nn and 2 H(e, e'7r _ )pp, which gave within the 
uncertainties a ratio of unity for the longitudinal cross 
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FIG. 2. Separated cross sections <tl and <tt (full and open 
symbols, resp.) compared to the Regge model (full curve 
for L, dashed curve for T). The Q 2 values are in units of 
(GeV/c) 2 . 

sections. Hence the VGL model is still considered to be 
a good starting point for determining F n . 

The comparison with the ctl data shows that the t de- 
pendence in the VGL model is less steep than that of the 
experimental data. As suggested by the analysis |l5j of 
older data, where a similar behaviour was observed, we 
attributed the discrepancy between the data and VGL to 
the presence of a negative background contribution to the 
longitudinal cross section, presumably again due to reso- 
nances. Since virtually nothing is known about the effect 
of these resonances on ctl, we proceeded on two paths to 
determine a trustworthy value of F n . First we fitted the 
VGL prediction for ctl to the data by adjusting F^ at the 
lowest \t\ bin, as shown in Fig. ^, where it is assumed to 
be most reliable, owing to the dominant t pole behaviour. 
However, since there is no reason to believe that the (neg- 
ative) background is zero at the lowest — t, the result is 
an underestimate for F n . Secondly, F n was determined 
adding a Q 2 dependent negative background to ax (VGL) 
and fitting it together with the value of F„. The back- 
ground term was taken to be independent of Z. This 
was suggested by looking at the 'missing background' in 
cr T , i.e., the difference between the data and VGL for 
ax. That background is almost constant or slightly ris- 
ing with \t\. Then, assuming that the background in ctl 
has a similar t-dependence, a constant background leads 
to an overestimate of F^. Our best estimate for F^ is 
taken as the average of the two results. The model un- 
certainty (in relative units) is taken to be the same for 
the four Q 2 points, and equal to one half of the average of 
the (relative) differences. The results are listed in Table 
[| and shown in the form of Q 2 F n in Fig. ^. The error 
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bars were propagated from the statistical and systematic 
uncertainties on the cross section data. The model un- 
certainty is displayed as the gray bar. The fact that the 
value of at Q 2 = 0.6 (GeV/c) 2 is close to the extrap- 
olation of the model independent data from j^] , and that 
the value of the background term is lower at higher W 
(see below) , gives some confidence in the procedure used 
to determine F v . 

For consistency we have re-analysed the older L/T 
separated data at Q 2 = 0.7 (GeV/c) 2 and W = 2.19 
GeV from DESY jll]]. We took the published cross sec- 
tions and treated them in the same way as ours. The 
background term in <tl was found to be smaller than in 
the Jefferson Lab data, presumably because of the larger 
value of W of the DESY data, and hence the model un- 
certainty is smaller, too. The resulting best value for F„, 
also shown in Fig. ||[ is larger by 12% than the origi- 
nal result, which was obtained by using the Born term 
model by Gutbrod and Kramer Here it should be 

mentioned that those authors used a phenomenological 
i-dependent function, whereas the Regge model by itself 
gives a good description of the f-dependence of the (un- 
separated) data from Ref. [jn^ . 

The data for F w in the region of Q 2 up to 1.6 (GeV/c) 2 
globally follow a monopole form obeying the pion charge 
radius ||. It should be mentioned that the older Bebek 
data in this region suggested lower F T values. However, 
as mentioned, they did not use L/T separated cross sec- 
tions, but took a prescription for ctt- Our measured data 
for <tt indicate that the values used were too high, so that 
the values for F n came out systematically low. 

In Fig. U the data are also compared to theoretical 
calculations. The model by Maris and Tandy |16| pro- 
vides a good description of the data. It is based on the 
Bethe-Salpeter equation with dressed quark and gluon 
propagators, and includes parameters that were deter- 
mined without the use of F^ data. The data are also 
well described by the QCD sum rule plus hard scattering 
estimate of Ref. M . Other models were fitted to the 
older F^ data and therefore underestimate the present 
data. Figure || also includes the results from perturba- 
tive QCD calculations ||. 

In summary, new accurate separated cross sections for 
the 1 H(e, e'7T + )n reaction have been determined in a kine- 
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FIG. 3. The Jefferson Laboratory and reanalyzed DESY 
values for fV in comparison to the results of several calcu- 
lations. The model uncertainty is represented by the gray 
area. The (model-independent) data from Ref. H are also 
shown. A monopole behaviour of the form factor obeying the 
measured charge radius is almost identical to the Maris and 
Tandy curve. 

matical region where the t-pole process is dominant. Val- 
ues for F^ were extracted from the longitudinal cross 
section using a recently developed Regge model. Since 
the model does not give a perfect description of the t- 
dependence of the data, our results for F^ contain a 
sizeable model uncertainty. Improvements in the theo- 
retical description of the 1 H(e, e'ir + )n reaction hopefully 
will reduce those. The data globally follow a monopole 
form obeying the pion charge radius, and are well above 
values predicted by pQCD calculations. 
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TABLE I. Best values for F„ from the present data and 
from the re-analyzed data from Ref. Jxij ] ■ The total (system- 
atic and statistical) experimental uncertainty is given first, 
and second the model uncertainty. 



Q 2 (GeV/c) 2 


W (GeV) 


Fir 


0.60 


1.95 


0.493 ± 0.022 ± 0.040 


0.75 


1.95 


0.407 ± 0.031 ± 0.036 


1.00 


1.95 


0.351 ± 0.018 ± 0.030 


1.60 


1.95 


0.251 ± 0.016 ± 0.021 


0.70 


2.19 


0.471 ± 0.032 ± 0.037 
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